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Abstract 

We analytically calculate the dominant two-loop electroweak correction, of 
0(G F mf), to the partial width of the decay of a Higgs boson, with mass Mh <C rut, 
into a bottom-quark pair, and describe the most important conceptual and tech- 
nical details of our calculation. As a by-product of our analysis, we also recover 
the 0{a s GFmD correction. Relative to the Born result, the 0(G F m^) correction 
turns out to be approximately +0.047% and, thus, more than compensates the 
0{a s GF'm't) one, which amounts to approximately —0.022%. 
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1 Introduction 



The standard model (SM) of elementary particle physics predicts the existence of a last 
undiscovered particle, the Higgs boson, whose mass Mh is a free parameter of the theory. 
The direct search for the Higgs boson at the CERN Large Electron-Positron Collider 
LEP 2 only led to a lower bound of Mh > 114 GeV at 95% confidence level pQ. On 
the other hand, high-precision measurements, especially at LEP and the SLAC Linear 
Collider SLC, were sensitive to the Higgs-boson mass via electroweak radiative corrections. 
These indirect measurements yielded the value Mh = (85^2s) GeV and an upper limit of 
Mh < 166 GeV at 95% confidence level [2]. The vacuum-stability and triviality bounds 
suggest that 130^M#^180 GeV if the SM is valid up to the grand-unification scale (for 
a review, see Ref. [3]). For these reasons, one hopes to discover the Higgs boson at the 
CERN Large Hadron Collider (LHC), which will be capable of producing particles with 
masses up to 1 TeV. The first question after discovering a new scalar particle will be if 
it actually is the Higgs boson of the SM, or possibly some particle of an extended Higgs 
sector. Therefore, it is necessary to know the SM predictions for the production and decay 
rates of the SM Higgs boson with high precision. Its decay into a bottom-quark pair is 
of special interest, as it is by far the dominant decay channel for M#<^140 GeV (see, for 
instance, Ref. [1]). 

At this point, we wish to summarise the current status of the calculations of radiative 
corrections to the H — > bb decay width in the so-called intermediate mass range, defined 
by Mw < Mh < 2M\y- The correction of order 0(a s ) was first calculated in Ref. [5]. 
The complete one-loop electroweak correction was found in Ref. [6]. As for the 0(a 2 s ) 
correction, the leading [7J and next-to- leading [8] terms of the expansion in m1/M H of 
the diagrams without top quarks are known. The diagrams containing a top quark can 
be divided into two classes. The diagrams containing gluon self-energy insertions were 
calculated exactly [9], while for the double-triangle contributions the four leading terms of 
the expansion in M H /m\ are known [10]. In Ref. [11], the 0{o? s ) correction without top- 
quark contributions was calculated in the massless limit. The correction induced by th top 
quark was subsequently found in Ref. [12] using an appropriate effective field theory. As 
for the correction of order 0(a s GFm 2 ), the universal part, which appears for any Higgs- 
boson decay to a fermion pair, was calculated in Ref. [13] and the non-universal one, using 
a low-energy theorem, in Ref. [H] . The latter result was independently found in Ref. [15] . 
Apart from the Higgs-boson decay into a it pair, only the one into a bb pair has such 
non-universal top-quark-induced contributions, as bottom is the weak-isospin partner of 
top. The universal and non-universal corrections of order Oi^Gp^t) were calculated in 
Refs. [16] and [17] , respectively. Finally, also a result for the universal correction of order 
0{G F mf) was published [TBI - 

In this paper, we calculate the complete correction of order O(Gpmf), including both 
the universal and non-universal contributions. To this end, we formally assume that 
Mh <C m t . This includes the intermediate mass range of the Higgs boson. Our result 
for the universal contribution in the on-mass-shell scheme agrees with the one found in 
Ref. [18J, after correcting an obvious mistake in the latter paper. The key results of our 
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calculation were already presented in a brief communication [19]. Here, the full details 
are exhibited. 

Our calculations are performed in 't Hooft-Feynman gauge. We adopt the on-mass- 
shell scheme and regularise the ultraviolet divergences by means of dimensional regular- 
isation, with D = 4 — 2e space-time dimensions and 't Hooft mass scale [i. We use the 
anti-commuting definition of 75. As a simplification, we take the Cabibbo-Kobayashi- 
Maskawa quark mixing matrix to be unity. The Feynman diagrams are generated and 
drawn using the program FeynArts [20] and evaluated using the program MATAD [2Tj . 
which is written in the programming language FORM [22]. 

In order to check our calculations, we also rederive the correction of order 0(a s GFm 2 ). 
Our result agrees with Refs. [T3"lll4|ll5] . Since this calculation follows the lines of the one 
leading to the 0(G 2 F mf) correction, being actually simpler, we refrain from going into 
details with it. 

This paper is organised as follows. In Section [2J we describe in detail the renormal- 
isation procedure underlying our analysis. In Section [31 we present the details of our 
diagrammatic calculations. In Section [U we explain how a part of our calculations can 
be checked through the application of a low-energy theorem. In Section [5], we evaluate 
the 0(G 2 F mf) corrections numerically and compare them with the 0(a s GFTn 2 ) ones. We 
conclude with a summary in Section [6j 

2 Renormalisation procedure 

For the reader's convenience, we present in this section the details of the renormalisation 
procedure which has to be carried out. We derive general expressions for the mass coun- 
terterms and wave-function renormalisation constants in the on-shell scheme, valid for 
any number of loops. Furthermore, we derive the tadpole renormalisation counterterms 
and describe the treatment of the corrections due to external legs. In our calculations, we 
do not need to consider electric-charge renormalisation constants, because, to the orders 
we consider here, there are no such contributions. 

Before going into details, we would like to mention that the expressions for the mass 
and wave-function renormalisation constants to be derived here are only valid for sta- 
ble particles. Instable particles do have complex self-energy amplitudes, so that their 
resummed propagators have complex poles. In that case, the renormalisation conditions 
are more complicated (see, for instance, Ref. [23J ) . Since all self-energy amplitudes ap- 
pearing in the calculations of this paper are real, we can restrict ourselves to the case of 
stable particles. 

2.1 Mass and wave-function renormalisation 

We write the bare masses in the Lagrangian as sums of the renormalised ones and the mass 
counterterms. In the on-shell scheme, we fix this splitting by the requirement that the 
renormalised masses are identical to the poles of the propagators including all radiative 
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corrections. Furthermore, the wave-function renormalisation constants are obtained as 
the residues of the propagators at their poles. 



2.1.1 Higgs-boson mass and wave-function renormalisation 

For the amputated one-particle-irreducible self-energy of the Higgs boson, we write 

" {l^l\- H --- = iX H (q 2 ). (1) 



Thus, the dressed propagator, including all radiative corrections, becomes 

^V) = + --{^---+—(^—(^—+... 



(2) 



q 2 -M 2 Hfi + Z H (q 2 y 

The on-shell renormalisation condition reads 

S H (M 2 H ) = 0. (3) 

Writing the bare mass of the Higgs boson as the sum of the renormalised mass and a 
counterterm, Mfj = Mfj + SM]j, we have 

5M 2 H = £tf(M|). (4) 

Here and in the following, it is understood that, in the expression for a counterterm, all 
bare quantities have to be replaced by the renormalised ones plus the respective counter- 
terms. In the case of the Higgs-boson mass counterterm, this means that E H (M|) has to 
be expressed in terms of renormalised quantities. For higher-order expressions, this has 
to be done iteratively. 

Expanding Eq. (T5]) about q 2 = Mfj and taking the limit q 2 — > Mjj, 

S H \q 2 )- ' ' 



q 2 — Mjj 1 + T! H (M 2 H ) + 0(q 2 - M 2 H 

> 



q 2 -M 2 H ' 

we read off the Higgs-boson wave-function renormalisation constant as 

7 1 
H 1 + E^(M|)- 



(5) 



(6) 
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Writing Zh = 1 + 5Zh and performing a loop expansion of Eq. ([6]), we have 

6Z$ = -£«'(il4), (7) 

5Z^ = -^\Ml) + {^\Ml)) 2 . (8) 

Here and in the following, numbers placed in parentheses as superscripts specify the loop 
order of the perturbative expression. 

2.1.2 Fermion mass and wave- function renormalisation 

The amputated one-particle-irreducible self-energy of fermion / has the form 



q 




l-PIJ * = iT, f (q) =i$u„E ftL (q ) + i$u + T< fjR (q ) + im/ l0 S/,s(5 ), (9) 



where 77i/,o is the bare mass of fermion / and a>± = (1 ± 7s)/2 are the projectors onto the 
helicity eigenstates. 

The fermion field / is composed of left- and right-handed components, I and r, re- 
spectively, as 

f = l + r, Z = w_/, r = w + /. (10) 

In the electroweak theory, I and r interact differently, which has to be accounted for in 
the renormalisation procedure. In terms of these components, the purely fermionic part 
of the SM Lagrangian reads: 

C = f{i$ — mffi)f = il$l + ir$r — mf^fl — mf t olr. (11) 

We see that I and r are massless fermion fields with propagators 

I r i 



q q 

In addition, we have the following r-l transition vertices: 

I r r I 



(12) 



irriffi. (13) 



From Eq. (jUJ), we read off the amputated one-particle- irreducible self-energies pertaining 
to the four different helicity combinations as 




i/l-PlV- = ^ fx (r- 



|^W)-^ = ^E /)J i(g a ), 

^J^y- = ~^Jl^y^- = im f>0 Z fjS (q 2 ). (14) 
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Note that above expressions do not yet include the tree-level contributions from Eqs. (jT2l) 
and (fT3|) . Equations (jT2l) - (TbfI) are the ingredients out of which we construct the propa- 
gators of the left- and right-handed fields including all radiative corrections. This is done 
in close analogy to the case of 7-Z-mixing (see, e.g., Ref. [24]). To this end, we introduce 
the propagator-type symbols 

I 




E Vpftti 



fj ^(l + £ /jL (g 2 ))' 




n=0 



E ^Mq 2 ) 



and the vertex-type symbols 




im f,0 ( S /,s(<? 2 ) - l) 




' 1 +^Yl-PlV- = ^ / , (S /i5 (g 2 )-l). 



(15) 



(16) 



Next, we evaluate the dressed propagator of the left-handed fermion field, including all 
radiative corrections, as 



I 



I ^ t ^ I I t _ I _ r ^ I 

^y^y^ + +^y+^y^y^y+ 



where 



^(l + S /iL (g2)) 
id 



l + S /iL (g2) f-m^ftf 



(17) 



;i - s /lS (9 2 )) 2 



(l + S / ,L(g 2 ))(l + S / , /? (g2))- v ' 

In a similar way, we find the dressed propagator of the right-handed fermion field, includ- 
ing all radiative corrections, to be 



s„ l (q) 



(19) 
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For completeness, we also resum the loop contributions by which a left-handed field con- 
verts into a right-handed one and vice versa. Proceeding similarly as in Eq. (TT7T) . we 
obtain 



<V(<?) 



I ^ r I ^ r ^ I ^ r 



OO r- 



X 



X 



im f;0 (l - S /iS (g 2 )) 



(20) 



" (l + E^tl + E^)) q 2 -mi f(q 2 y 

Since Eq. ( 1201) is symmetric under the interchange of the indices L and R, we also have 

Stf(q) = S?(q). (21) 

We now derive the fermion mass counterterm. Writing mf t o = mf + 5m/, where mj 
is the renormalised mass and 8m f is the mass counterterm, and imposing the on-shell 
renormalisation condition, 



S ij (q)u f (q)\ c 



0. 



where ij = 11, rr, It, rl and Uf(q) is the spinor of the incoming fermion /, we obtain 

5m f 1 



m f 



f K) 



I. 



(22) 



(23) 



Expanding Eq. 



5m^ 



we find the explicit one- and two-loop expressions, 



^Umi) + mUm 



m f 
Sm^p 
m f 



(i) 



,(2) 



mi) + ^Umi) + £ v f %(m 



f,S\ nb f)i 



(2) 



(24) 



5m 



(i) 



(25) 



The one-loop expression of Eq. ff24l is well known (see, e.g., Ref. [6]). The two- loop 
expression of Eq. ( l25l) agrees with the one obtained in Ref. [25] using an alternative 
procedure. 
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Finally, we derive the wave-function renormalisation constants for the left-handed and 
right-handed fields. Expanding Eqs. (ITT)) and (TT9T) (|2TT) about d = mj and taking the 
limit d —>■ rrif, we have 



S ll/rM) 



1 



q 2 — m 2 



f (l + Z f>L/R (m} 



-i 2 /'( m /) 

/ f(mj) 



+ 0(q 2 - m 2 ) 



irrif 



m 



2 ' 

f 



q 2 ^m) im f ^/Z f)L Z fA 



2 f'(mj) 



m. 



+ 0(q 2 - mj) 



(26) 



where 



J f,L/R 



(1 + Z f , L/R (™})) (l - m}^ 



(27) 



Writing Zj L / R = 1 + 5Zf^ L / R and performing a loop expansion of Eq. (|2T)) . we have 





y(l) 


y(l)' 
~~ L 


- Eg" 


- 2Eg", 




-Eg) 


y(l)' 
_ L 


- Eg" 


oy(l)' 




y(2) 
L 


y(2)' 
_ L 


- Eg" 


- 2Ef + E« 



l +24 1); + S R 

2 



+ S «4 1)/ -24 1) 4 1)/ + (4 1)/ + 4 1)/ + 2s 



+ 2£f 
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(2) 

f,R 



,(2) 



y (2)/ y (2)/ 9V (2)/ 
L ^R ZZj 5 



+ e« (4 1) + 4 1)/ + 24 1)/ + 24 1) ' 



+ 4 1) 4 1)/ -24 1) 4 1)/ 



2S 



(i)< 



Here, we used the abbreviations 

s (n) 



(n)/ 
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(28) 
(29) 

(30) 
(31) 



(32) 



where X = L,R, S. These expressions again agree with Refs. 

If parity was conserved, we would have 4(g 2 ) = ^i?(^ 2 ) anc ^ thus recover the structure 



iZf 



m f 



(33) 



which is familiar from quantum electrodynamics. 
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2.1.3 W-boson mass renormalisation 



The amputated one-particle-irreducible self-energy of the W boson can be decomposed 
into a transverse and a longitudinal part as 




(34) 



where 



A"" = g 



[IV 



q H q u 



q 



q»q v 
q 2 



(35) 



Owing to the loop-induced mixing of the W boson with the charged Higgs-Kibble ghost (p, 
we must also take into account the one-particle-irreducible W <-> <fi transition amplitudes 
and the one-particle-irreducible 0-boson self-energy, 



l-PlV = iq^w^q 2 




-iq^ W(P (q 2 ), 



^-{^^---^^(g 2 ). 



(36) 



In 't Hooft-Feynman gauge, the bare propagators of the W and bosons are given by 



Mb* 



n 2 _ 2 ' 

q m w,o 



(37) 
(38) 



with a common bare mass M\y )0 . In order to obtain the dressed IV-boson propagator, we 
proceed in two steps. In the first step, we resum the one-particle irreducible self-energies 
of the W and bosons separately. In the second step, we systematically combine these 
results by accommodating all possible W <-» <fi transitions. 

The resummation of the one-particle irreducible P^-boson self-energy leads to 



W 



w 




+ ^4 i_pi KX/v + Xkl i_pi KXi i_pi KX + 



G w ,^{q 2 ) + G w ^ a (q 2 ) (-ill$(q) \ G w ^ u (q 2 ) 

+ G w ^ a (q 2 ) (-in^(g)) G w ^(q 2 ) (-ill^(g)) G WM {q 2 ) + .... (39) 
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The series in Eq. (|39|) may be resummed by inserting Eqs. (1341) and (|37|) and exploiting 
the identities 




Gw,fia(<l 



A"„A" P = A%, 
A^q vp = 0, 



A a u Z w>T (q 2 )+q a ^ w>L (q 2 



(40) 



g 2 -M 2 



L w,o 

A\(E w , T (q 2 )) 2 + q\^wAq 2 )f 



+ 



(q 2 - Ml y 



Gw,fj.a(q 2 



G 



g\ + A\ 

g\ + A a 



I ^2 _ /\f 2 

n=1 \ 9 M w,o 
1 



1 S w , T (q2) 

1 ' a 2 -M 2 

1 m W,0 



! I + q a 



-z w , L (q 2 ) 

i 



Sw-,L(g 2 ) 

9 M W,0 



- 1 



1 + „u_ w" 



A 



q 2 - M*r )0 + ^wAq 2 ) q 2 - m^ )0 + Z w Aq 2 ) 

(^w!pure) ^ (<?)• 



(41) 



The resummation of the one-particle-irreducible 0-boson self-energy proceeds in analogy 
to the Higgs-boson case discussed in Section 12.1.11 and yields 



-'-0 



- - -f i-pi V - -f i-Pi V --+... 



g 2 - M*, t0 + E^g 2 



(42) 



The contribution of unmixed PF-boson propagation in Eq. (1411) needs to be comple- 
mented by the contribution that emerges by combining it with the contribution of un- 
mixed 0-boson propagation of Eq. fj42|) via the one-particle-irreducible W <-> <p transition 
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amplitudes in all possible ways. This additional contribution is given by 

, , , . . W A ± A W w^±^w^t^w 
(<V,mixj Mi , (?) = — -©-x^- + -^4^ _ ^P>™^ -_- ^Pjh^ + . . . 



q^w4>(q 



2\ 



g 2 - M^ Q + S w ,L(g 2 ) q 2 ~ M^ + £ (g 2 ) 

g 2 (W? 2 )) 2 



oo 

X 

n=0 



X 



(g 2 - M^ + £^(g 2 )) (g 2 - M^ + £ (g 2 )) 
-g„£ w<)! ,(g 2 ) 



g 2 - + Z w , L (q 2 ^ 

~ iq >" ~ q2 (43) 



g 2 - M^ + S^ iL (g 2 ) , _ (g 2 -M^ o +E w , L (^))( g 2_ M ^ o+S0((?2 )) 

Adding Eqs. (j4"Tl) and fj43]) . we obtain the fully dressed VT-boson propagator as 

( S w), u (q) = + (^mix)„ (?)• (44) 

Its inverse is found to be 

= ^(? 2 - ^, ) + *A^ T (g 2 ) + W Uw^q 2 ) - /f/^] 2 \ ■ (45) 

\ q ~ m w,o + 2j </>w j j 

The on-shell renormalisation condition reads 

S^(q 2 )e w , u (q)\ q2=Mh =0, (46) 

where e^(g) is the polarisation four- vector of an external W boson. Writing M^ = 
+ and exploiting the transversality property g M e^^(g) = 0, we finally have 

5M^ = Z w , T (M%r). (47) 

We note in passing that Eq. (j4"Tj) is not influenced by W <-> mixing. 



2.2 External-leg corrections 

In this section, we discuss the structure of the amputated matrix element A for the decay 
process H — > bb and explain how to obtain from it the transition matrix element T by 
incorporating the wave-function renormalisation constants. 
The general form of A reads 




= i (A + i x A 2 + f 2 Aa + i 2 i x A* + 7sA + 7s^ A + 7s0 2 A + 75^ A) , (48) 
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where q\ and q 2 are the four-momenta of the outgoing b and b quarks, respectively, and 
Ai (i = 1, . . • , 8) are scalar form factors. Projecting onto each of these form factors, we 
observe that, to the orders we consider in this paper, only two of them are independent. 
In fact, we have 

A 2 = -Aa = A 6 = -At, 

A 4 = A 5 = A 8 = 0, (49) 
so that A collapses to the simple form 

A = A A + A B (j 2 - i^j to- , (50) 

where A a = Ai and Ab = — 2A 2 . 

Then, T is obtained by dressing A with the renormalised wave functions of the external 
legs as 

T = ^fz~^ Z btR u r (q 2 , r 2 ) + ^/Z bjL ui(q 2 , r 2 )J A (^/ Z b)R v r (qi, n) + ^/Z b ~J / Vl(q 1 ,r 1 )^ 
= y^u b (q 2 , r 2 ) {^JZ h ^ R u]- + a/ Z biL Ld + ^j A \ J Z biR u + + Z b)L uj^j v b {q x , n), (51) 

where v b {q\,r\) and u b (q 2 , r 2 ) denote the spinors of the outgoing b and b quarks with spins 
ri and r 2 , respectively. Inserting Eq. fl50|) into Eq. floTT) . we obtain the master formula 

T = \[Z~^ \y/ Z biL Z biR A A + m b Z btL A B ^ u b (q 2 , r 2 )v b (q u ri). (52) 
Note, that the terms involving 75 vanish upon application of the Dirac equation. 



2.3 Tadpole renormalisation 

As is well known (see, for instance, Ref. [26]), one can introduce a so-called tadpole 
renormalisation in order to avoid the calculation of diagrams containing tadpoles. For 
the reader's convenience, in this section, we rederive the counterterm vertices of the 
tadpole renormalisation along with the counterterm vertices of the Higgs-boson mass 
renormalisation. 

The tadpole renormalisation concerns only the Higgs part of the SM Lagrangian, 

£ H i g gs = (D^y(D»$) + fi 2 &$ - ^($ f $) 2 , (53) 

where $ is a weak-isospin doublet of two complex scalar fields. The free parameters, \i 
and A, are chosen in such a way that one stays with a non-vanishing vacuum expectation 
value v, which is defined by 

^ = |(0|$(x)|0)| 2 = ^. (54) 
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If we parameterise 



and substitute fi and A by 



<I>( ^ i i / t + /f • / ^ ) ( 55 ) 



2 



4 



= V + (56) 



Eq. (1531) takes the form 



1 1 M 2 

Cm ggs = -(D»H)(D»H) + ~(D» X )(D» X ) + (D^D^) + tH- -fH 2 



x 2 



(57) 



where = (</> + )^. We see that Mh has the physical meaning of the Higgs-boson mass. 
In this step, we did not exploit Eq. fl54l . which implies that t = 0, so that we could just 
have emitted all terms containing t. However, as was argued above, it is useful to keep 
them and to renormalise t along with M H by substituting 

t -> t = + St, 
M 2 H -> M 2 Hfl = M 2 H + 5M 2 H (58) 

in Eq. (|5"T|) . Notice that Eq. (|57|) represents a bare Lagrangian, so that v, t, and Mh are 
actually bare parameters. For consistency, we thus also substitute v — > Vq. Then, Eq. (J5T 
becomes 

1 1 M 2 

£m gg s = -(D,H)(D»H) + -{D,x){D»x) + (IV) W") " -y^ 2 

- (tf 2 + x 2 + 20-0+) - ^ + X 2 + 20-0 + ) 2 

2 2w 2v V^o / 

x (H 2 + X 2 + 20-0 + ) " A f- + (H 2 + X 2 + 2r + ) 2 



(59) 



From the terms proportional to St and SM H , we can read off the desired counterterm 
vertices, which we list in Table HJ 
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Table 1: Counterterm vertices related to the Higgs-boson tadpole and mass renormalisa- 
tion. 



H: 



HH: 



XX- 



iSt 



-i8M% 



.St 
i — 

v 

St 

5 — 



HHH: 

HXX- 

H(/><j>: 



3 

-i — 
1 

-% — 

v 



St e „ 

h SM 

vo 



- + 5M 2 
vo 



St 



-i— — + SM Z H 

V \Vo 



HHHH:- 

xxxx- - 

HHxX- - 
HH(f)(j): - 

xxH- - 



-l-a — 



-l-a — 



+ 8 M 2 H 



+ 8M 2 H 



SM 2 H 



+ SM 2 H 



3 l- + SM H 



+ SM 2 H 



The Higgs-boson mass renormalisation condition was already discussed in Section |2TTTTJ 
As a renormalisation condition for St, we set 

St = — T, (60) 
where T stands for the sum of all amputated one-particle-irreducible tadpole diagrams, 




As can be seen from Table [TJ there is a one-point Higgs-boson counterterm vertex, iSt, 
that forces a cancellation with all diagrams having a tadpole at its place. Therefore, 
upon tadpole renormalisation, one does not have to consider tadpole diagrams anymore. 
However, now one has to take into account all the tadpole counterterm vertices in Tabled], 
except for the one mentioned above. 



3 Results 

In this section, we present the details of our actual calculations. After making some 
general remarks, we describe in Sections 13.11 I3.2[ and 13.31 the explicit computation of 
the decay rate at tree level, at the one-loop order 0(Gftti 2 ), and at the two-loop order 
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0(G F mf), respectively. Section [37T1 also contains the expressions for the renormalisation 
constants at order 0(Gpm^), which are needed in the one-loop and two- loop calculations. 

In order to compute the leading large-m 4 contributions of the various two-loop dia- 
grams, we apply the asymptotic-expansion technique (for a careful introduction, see 
Ref. |27J). However, it turns out that all non-trivial contributions of the self-energy 
and Hbb vertex diagrams (see Figs. HI O and [6]), which are of leading order in mt, cancel 
among themselves or, in case of the W^-boson self-energy, in combination with complete 
counterterm diagrams arising form the Higgs-boson tadpole and mass renormalisations. 
Specifically, in Fig. HJ there are non-naive contributions due to the asymptotic expansion 
of diagrams (i)-(o) that cancel against diagrams (p)-(v); in Fig. the non-naive contri- 
butions of diagrams (a) and (t) cancel; and in Fig. [6] those of the diagrams (e) and (i) 
cancel. After these cancellations, only naive contributions due to diagrams involving top- 
quark propagators remain. Therefore, we can naively expand in all masses and momenta 
except for the top-quark mass and retain only the leading terms. Obviously, this requires 
the Higgs-boson mass to be smaller than the top-quark mass, which is compatible with 
the intermediate-mass range of the Higgs boson, as mentioned in the Introduction. 

The ultraviolet divergences which have to disappear in the final expression for the 
decay rate are cancelled through the application of the renormalisation procedure, which 
we carry out in the on-mass-shell renormalisation scheme. This provides a non-trivial 
check for our calculations. As explained in Section 12.3} we use the counterterm vertices 
of Table [T] for the Higgs-boson tadpole and mass renormalisations. However, while we 
renormalise the Higgs-boson mass already at the Lagrangian level, we replace all other bare 
parameters at the end of the calculations without recourse to any counterterm vertices. 
This procedure turns out to be most convenient for our purposes. 

As a further check on our calculations, we also rederive the correction of order 
0(a s Gpm1). This result is presented in Section 13.41 Finally, we apply a Higgs-boson 
low-energy theorem [28J, which allows for an independent calculation of the various Hbb 
diagrams at order 0(G F mf). This is explained in Section HI 



3.1 Tree-level result and O(GFm^) renormalisation constants 

The tree-level diagram is depicted in Fig.[T](a). Using the notation introduced in Eq. ( 1501) . 
the corresponding amputated matrix element is in bare form written as 

(62) 

The tree-level transition matrix element is 

T(°)=4 0) n fe (g 2 ,r 2 )^(g 1 ,r 1 ), (63) 

and the decay rate is 

(0) = V2N c G F M H ml f _ 4mf \ 3/2 

Svr V Ml) ' 1 J 
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Figure 1: Diagrams contributing to H — > bb at (a) tree level and (b) order ©(Cr^raf). 

where N c = 3 is the number of quark colours. Furthermore, we have introduced Fermi's 
constant Gf via the Born relation 

- = 2 1 / i G 1 J 2 . (65) 

v 

In the following, we have to renormalise the vacuum expectation value. Through the 
order of our calculations, this can be achieved by writing [29] 

- = 2^Gf Q , (66) 



v 

with 



M 2 W 



G F ,o - G F —-^—. (67) 

m w,o 

Thus, the renormalisation of the vacuum expectation value is reduced to the one of the 
jy-boson mass. 

In the remainder of this subsection, we list all relevant renormalisation constants of 
order QiGprn^). They are derived by evaluating the diagrams of Fig. [2] and applying 
Eqs. flU), dZD, (JMD, (HI, (EU), (JHD, and (JS0D- Since we shall compute the correction of 
order 0(G F m^), these renormalisation constants are needed through order 0(e) in the 
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expansion in e. The results read 

(Jfto = C efi x tfi ml v N c 



5Z 



5m 



(i) 

H 

(1) 



5Z, 



(i) 

6,L 



5m 



(i) 



m t 



C tfi x tfi ml Q N t 
0, 

C £!O x tfi Ml N ( 



- + 4+(4 + 2C(2))e + 0(e 2 

.^_ 4 + (_4_ 6C (2)) e + 0( e 2 
e 



-- + |-C(2)e + (9(^ 
e 3 



3 5 
27 ~ 4 
1 3 
7 ~ 2 " 



-|-? C (2))e + 0(e* 
•^-^(2))e + 0(e 2 ) 



| + 4+(9-^C(2))e + 0(^ 



2 

1 

e 



---C(2))e + 0(^ 



where we use the abbreviations 



a 



/ 4tt^ 2 ( 

V m t ( 

G F m\ 
8tt 2 V^' 



-7b 



(68) 
(69) 
(70) 

(71) 

(72) 
(73) 
(74) 

(75) 



(76) 



with 7£ being Euler's constant. 



3.2 Correction of order 0{Gpf^) 

At order 0(Gprnj), only the one diagram depicted in Fig. QJb) contributes. Using the 
notation of Eq. (1501) . we obtain for the expansion in e through order 0(e): 



"^4 ~~ < ~ y e,02^,0 

^0 



- + 2 + (2-C(2))e + (9(6 2 ) 
e 



A%^C e>0 x t>0 - f-l-|e + 0(e 2 ; 

Vq \ Z 



(77) 



Expanding Eq. fl52l) and replacing the bare masses by the renormalised ones plus their 
counterterms in Eq. f[6"2"j) . we find the transition matrix element to be 



(1) 



5m u 



(i) 



(78 
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Figure 2: One-loop self-energy and tadpole diagrams contributing at order 0[GFW? ti 



where is the amputated matrix element of Eq. (|62|) and 



K ~- Zh 2 Ml 



2(1) 
W 



(79) 



w 



is the one- loop contribution to the universal counterterm 5 U , which exhausts the full 
OiGpW^) corrections for Higgs-boson decays to fermion-antifermion pairs, except for 
those into ti and bb pairs. For simplicity, we omitted the spinors on the right-hand side 
of Eq. (1781) ; we shall also do this in the following. 5u~* and are ultraviolet finite and 
read 



# } = XtN t 



7 



7™ = T^xt [ -3 + N c 7 - 



The C(Gi?m^) correction to the decay rate thus becomes 



(80) 
(81) 



x+ — 6 + N c - 
p(o) M c 3 



(82) 



where is given in Eq. fl64l) . The results of this subsection are in accordance with 
Ref. ©. 
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3.3 Correction of order O(G^ra^) 

Expanding Eq. (1521) up to the two-loop order and replacing all bare masses in the tree-level 
and one-loop amputated matrix elements by the renormalised masses plus the correspond- 
ing counterterms, we find the following master formula for the transition matrix element 



(2) 



.(1) 



5m 



(i) 



(i) 

B,0 



^ (0) 



+ 



1 5m 



(i) 



2 m b 



m b 



+ *v 



^+2(1 



5m 



(i) 



m t 



-5Z {2) 
2 b ,R " 



SZ. 



(i) 

6,£ 




where 



2(2) 
W 



M 2 W 



~ (*& 2 



2(1) 



W 



3 / 5M, 



2(1) 
W 



M 2 W 



is the universal counterterm. 



3.3.1 Universal counterterm 



Let us first calculate the universal counterterm. To this end, we need the two-loop ex- 
pressions for SZh and SM^. The unrenormalised expressions are obtained by evaluating 
the diagrams in Figs. [3] and H] and applying Eqs. (jSJ and (jlTj) . the results being 



5Z 



(2) 
H,0 



0(e) 



= C>t 2 o^ iV c (i + i-f + 17C(2) + 0(e)) 



55) 



in accordance with Ref. [18J. In addition, there are contributions from the renormalisations 
of the bare parameters in Eqs. (1701) and (175]) . so that 



(2) 
H 



2(1 



5m 



(i) 



5M 



2(1) 



m t 



M 2 W 



5Mf ] = 6M$> + 2(1 



L w,o 



(86) 
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Figure 3: Higgs-boson self-energy diagrams contributing at order 0{G\ 



ml 



We are now in a position to specify the universal counterterm at order 0(G 2 F mf) as 
defined in Eq. ( l84|) . The result is 



*?=^c(y-6C(2) + ^) 



x+ 



495 



3tt< 



*7) 



If we convert Eq. flHTj) to a mixed renormalisation scheme which uses on-shell definitions 
for the particle masses and the definitions of the modified minimal-subtraction (MS) 
scheme for all other basic parameters, then we find agreement with Eq. (15) for x = in 
the paper by Djouadi et al. [IS] . However, the corresponding result for the pure on-shell 
scheme presented in their Eq. (27) for x = disagrees with our Eq. (j57|) . We can trace 
this discrepancy to the absence in their Eq. (25) of the additional finite term 5u^ Ap^ 
which arises from the renormalisation of the one-loop result in their Eq. (7) according to 
the prescription in their Eq. (18). 



3.3.2 Complete transition matrix element 

(2) 

Having provided o« , we now turn to the residual ingredients entering the transition 
matrix element of Eq. fl83l) . Evaluating the Hbb diagrams shown in Fig. [5j we find the 
form factors in Eq. (150]) at order 0(G^m^) to be 



A 



(2) 
A,0 



V 



» 2 r 



c,0 



±-5-5 

e l e 



7C(2) + N c ( ~ - - - 14 - 2C(2) 1 + 0{e) 



A 2 ) _ n 

v 



2 

e,0 



- + 1 + N c - + 9 



0(e) 
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Figure 4: VT-boson self-energy diagrams contributing at order 0(G 2 F mf). Insertions of 
—iSMjj in Higgs-boson lines and of i5t/v in 0- or x-boson lines are indicated by crosses. 
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Figure 5: Diagrams contributing to H — > bb at order 0(G 2 F mf). Insertions of iSt/vo in 
0-boson lines and of —i (5t/v + SMjj) /v Q in Hcpcp vertices are indicated by crosses. 



Evaluating the diagrams depicted in Fig. and using Eqs. f[2"5"j) . f[3"Uj) . and (jHTl) . 
we obtain the bottom-quark mass and wave-function renormalisation constants at order 
0(G 2 F mf). The renormalisation constants in terms of bare parameters read 



6m ?l 

m b 


- r 2 r 2 


w(2) 
b,L,0 


- r 2 r 2 


OZj b,R,0 


= 0. 



'27 31 13 59 _ A / 3 15 55 3 
2 7 l + 6C(2) + iV c fI + l + ^-C(2)Vw 



0(e) 



2e 



2e 



39) 



Additional contributions arise from the replacement of the bare t-quark and FT-boson 
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Figure 6: 5-quark self-energy diagrams contributing at order 0(G F mf). Insertions of 
i5t/vo in 0-boson lines are indicated by crosses. 



masses in Eqs. (17TT) . ( 1721) . and ( 1731) . so that 



Sm {2) 5m[ 2 ,j 5 m 



+ 

m b m b m b 

0/j b,L - OZl b,L,0 + b,L 

6zi% = 0. 



2(1 -el 



5m 



(i) 



6M, 



2(1) 



2(1 



I — — 

m t 



M 2 W 



6M, 



2(1) 



M 2 W 



(90) 



Now all ingredients for the evaluation of the renormalised transition matrix element 
of order 0(G 2 F m^) according to Eq. ( 1831) are available. We find 



T (2) = T (0) x 2 



2Q 49 
-y + iV c (18-6C(2)) + iV c 2 - 



(91) 



Adding Eqs. ( 1631) . ( IHTT) . and ( 19TI) . squaring, and extracting the 0{G\raf) term, we 
have 



r(2) 



r (o) •** 

= x 2 (116-6tt 2 ). 
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20 + iV c (29- 12C(2)) + AT C 2 — 

9 



(92) 



3.4 Correction of order 0(a s GFm^) 

As a by-product of our analysis, we can also compute the 0(a s G F m 2 ) correction to the 
H — > bb decay width. The comparison of our result with the literature [TBI 141)1 5] provides 
a partial check of our 0(G 2 F m^) results. Note, however, that the calculation considerably 
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simplifies as one passes from order 0(Gpfnf) to order 0(a s G F m 2 ). Using our tools, we 
indeed recover the well-known 0(a s G F mf) results for the universal correction [13] and 
the correction to the H — > bb decay width ^4]fl5] . 



x(x t a s ) 

p(X t Q! s ) 

r(°) 



X t ^C F N r _ 



71 



7T 



3 
4 



C(2) 



-12 



91n^ + iV 

M? + c 



15 

T 



C(2) 



7 Ml 
2 M, 2 



(93) 



respectively, where X t = G F M?/ (8n 2 V2) and C F = (iV 2 - l)/(2N e ). In Eq. (j23J, the 
bottom- and top-quark masses are denoted with capital letters, Mb and M t , respectively, 
to indicate that they are pure on-shell masses, i.e. they are defined in the on-shell scheme 
also with regard to quantum chromodynamics (QCD). The obvious disadvantage of this 
choice is the appearance of large logarithms of the type In (M^/m 2 ) starting already in 
order 0(a s ), which spoil the convergence behaviour of the perturbation expansion. As is 
well known [5], these logarithms can be resummed into the running bottom-quark mass, 
if m;, appearing in Eq. (1M1) is QCD-renormalised in the MS scheme at scale // = Mh, by 
substituting m b = m&(Mff). For consistency with the 0(G F m 2 ) and 0(G%mf) results 
presented above, which all refer to the electroweak on-shell scheme, we continue our 
discussion in a mixed renormalisation scheme where the on-shell definition of bottom- 
quark mass is adopted for electroweak corrections and the MS one for QCD corrections. 
Since we wish to treat the masses of the top and bottom quarks on the same footing, 
we adopt this mixed scheme for the top-quark mass as well. Furthermore, the analysis 
at order (9(a 2 GVm 2 ) [TBIfTT] reveals that Eq. f[9"3"j) is further improved according to the 
renormalisation group if m t and a s are taken to be m t = m t (m t ) and a s = ai n/ ^(m 4 ) with 
rif — 6 quark flavours, respectively. In this improved renormalisation scheme, Eq. fl93l) 
takes the form 



r(°) 



x t -C F N c 

7T 



19 
12 

a s /19 7T 2 

Xt ~ It ~ t 



C(2) 



x t — C F 

7T 



/I 57 

- 36 + A Ht7" C(2) 



o. 



7i 



13 

T 



-7T 



(94) 



To the order considered here, we have 

m t = M t \l- 



af\M t ) 



Cf. 



7T 



(95) 
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4 Low-energy theorem 



In this section, we present an alternative way of calculating all but one of the Hbb diagrams 
at order 0(G 2 F mj) which is based on the Higgs-boson low-energy theorem [28J. In fact, 
the Hbb diagrams of Fig. [5j with the exception of diagram (t), can be generated from the 
bottom-quark self-energy diagrams of Fig. [6] by in turn attaching an external Higgs-boson 
line to each of the top-quark lines. Diagrammatically, this can be represented as follows: 

\H 

t(q) — > t(q) ! t(q) 



imt >° ' (96) 



Here, we also made use of the fact that, in the large- m t approximation, the external Higgs 
boson does not carry any four-momentum into the respective diagram. Thanks to the 
identity 

% -im tj0 i m t ,o d f i \ 



the amputated matrix element of H — > bb is in the large-m^ limit related to the bottom- 
quark self-energy as 

Ao = ^^-^ (98) 

where it is understood that the differential operator only acts on masses which stem from 
propagators, not to those occurring in vertices, and that all quantities in Eq. (1981) are 
taken to be bare. Exploiting the structures underlying Eqs. (Q and fl50|) . Eq. fl98|) can be 
decomposed into two scalar equations. Identifying the four-momentum q in Eq. (Q with 
q 2 in Eq. fl5Dl) and noticing that q 2 = —qi in the soft-Higgs limit, we have 

A - m *>° 9 V 
•AAfl — m b,0 7^ ^b,Si 

A B '^ » Tf*. (99) 
2 v dm tfl 

The fact that the H — ► bb amplitude does not contain a term proportional to — ^ 1 )uj + 
is reflected by the fact that the right-handed part of the bottom-quark self-energy, 
vanishes to the orders considered in this paper. 

The results for Aa,o and Ab,o obtained through Eq. (I9~9l indeed agree with the direct 
evaluation of the respective diagrams in Fig. [5j 

5 Numerical results 

Finally, we explore the phenomenological implications of our results. Adopting from 
Ref. |S03 the values G F = 1.16637 x KT 5 GeV" 2 , af\M z ) = 0.1176, M z = 91.1876 GeV, 
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Table 2: Numerical values of the relative corrections A| , , A q , and r^/T^ to the 
H — > l + l~ , H — > qq, and H — > bb decay widths, respectively, at orders x = Gpm 2 , G 2 F m\, 
and a s G F m 2 . 



Order x 


A\ x) A { q x) rw /r(°) 


0(G F m 2 ) 
0(G\m\) 
0(a s G F m 2 ) 


+2.021% +2.021% +0.289% 
+0.064% +0.064% +0.047% 
+0.060% +0.452% -0.022% 



and M t = 174.2 GeV for our input parameters, so that a s (m t ) = 0.1076 and m t = 
166.2 GeV, we evaluate the relative corrections /T^ to the H — > bb decay width to 
orders x = G F m 2 , G 2 F mf, and a s Gpm 2 . For comparison, we also evaluate the relative 
corrections to the H — > l + l~ and H — > qq decay widths, where I = e, //, r and q = u,d, s, c, 
which, to the orders considered here, are given by 

a, = (i + <y 2 - 1 

= 2^)+24 2) + (4 1) ) 2 + 2e Qs) , 
A 9 = (l + A QCD )(l + <y 2 -l 

= A QCD + 25^ + 25^ + (5«) 2 + 25^ + 2A QCD 5i 1) , (100) 

where [5] 

a 17 

Aqcd = -C F -^ (101) 

7T 4 

is the 0(a s ) correction in the limit m q <C Mjj, with m q = m g (M^). 

The results are listed in Table [2j We observe that the 0(G 2 F mf) correction to 
increases the enhancement due to the 0(Gpm 2 ) one by about 16% and has more than 
twice the magnitude of the negative 0(a s Gpm 2 ) one. 

6 Conclusions 

We analytically calculated the dominant electroweak two-loop correction, of order 
0(Gpmf), to the H — >• bb decay width of an intermediate-mass Higgs boson, with 
M H < m t . 

We performed various checks for our analysis. The ultraviolet divergences cancelled 
through genuine two-loop renormalisation. Our final result is devoid of infrared diver- 
gences related to infinitesimal scalar-boson masses. We reproduced those Hbb triangle 
diagrams where the external Higgs boson is coupled to an internal top-quark line, which 
we had computed directly, through application of a low-energy theorem. After switching 
to a hybrid renormalisation scheme, our 0(G 2 F mj) result for the universal correction S u 
agrees with Ref. [18 J . Using our techniques, we also recovered the 0(a s Gpm 2 ) correction 
to the H — > bb decay width as well as the universal correction 5 U in this order. 
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The 0(G 2 F mf) correction to the H — > bb decay width amplifies the familiar enhance- 
ment due to the OiGpm^) correction by about +16% and thus more than compensates 
the screening by about —8% through QCD effects of order 0(a s GFmf). 
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